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Abstract 

We use an optimised perturbation expansion called the linear <5-expansion to 



> 

, study the phase transition in a Higgs sector with a continuous symmetry and large 

couplings. Our results show how to use this non-perturbative method successfully 
for such problems. We also show how to simplify the method without losing any 
, flexibility. 

On 

Quantum Field Theory has had many successes over the last fifty years but most 
are based on the use of small coupling perturbation expansions. However, there are 
Ph! many interesting problems where such an expansion is not appropriate. QCD is the 

classic example but phase transitions at non-zero temperatures even in models with small 
coupling constants, e.g. the Electro Weak model, are now known to be in this class. It is 



this latter example which motivates our work which will focus on the non-perturbative 
behaviour of the Higgs sector. 

Alternatives to perturbation theory, non-perturbative methods, have many limita- 
tions. For instance it is usually very difficult to go beyond the lowest order in methods 
such as large N [1] while numerical Monte Carlo codes are relatively expensive [2, 3]. An- 
other method with a long history is the Linked Cluster Expansion [4, 5] which now also 
exploits large amounts of computing power [6]. In this letter we extend a non-perturbative 
method called an OPE (Optimised Perturbation Expansion). In this approach one ex- 
pands around a solvable model which contains several arbitrary parameters, and thus is 
a variational ansatz. The expansion is not necessarily in terms of some small parameter. 
One then chooses the variational parameters, invariably using some optimisation criteria 
which is a highly non-linear procedure. The optimisation criteria has to ensure good 
results from a few terms of the series but there has been much discussion about what is 
a good optimisation criteria. 

As described, the OPE is a very general method so not surprisingly it has been re- 
discovered several times in different guises, applied to very different applications, and 
appears under many different names; it is also known as the linear ^-expansion [7], 
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action- variational approach [8], improved gaussian approximation [9], variational per- 
turbation theory [10], method of self-similar approximation [11], or the variational cu- 
mulant expansion [12]. For instance the method has been to the evaluation of simple 
integrals [13, 14, 7], solving non-linear differential equations [15], quantum mechanics 
[16, 17, 18, 9, 7] to quantum field theory both in the continuum [9, 10] and on a lattice 
[7, 19, 20, 21, 22, 23, 24, 25, 26, 27, 8, 12]. 

To motivate this work we first note that the analysis of pure gauge theories on the 
lattice using OPE used a variety of optimisation schemes when fixing the variational 
parameters. One can minimise the free energy [22] or demand that the series converges 
as fast as possible by minimising the high order terms with respect to the lower order 
terms [8] - the principle of fastest convergence. However, we believe the best results for 
pure gauge models are obtained when using the principle of minimal sensitivity [28], as 
discussed in [7, 19, 20, 21]. In this case the variational parameters, say {a}, are set 
separately for each physical quantity, say O, under consideration by demanding that the 
quantity changes as little as possible if the variational parameters {a} are varied from 
their optimal values {a}, specifically 



d(0) 



daj 



0. (1) 



Combinations of these criteria are also possible [11]. 

Given the interest in non-perturbative effects in gauged Higgs models it is interesting 
to note that the situation for the OPE in these cases [23, 24, 25, 26, 27] is not nearly 
so good, and the results are not as extensive as for the pure gauge case. The analysis 
of the gauged U(l) Higgs model [23, 24] fixed the variational parameters by minimising 
the lowest order calculation of the free energy, a procedure which is justified by noting 
the identity (exp{— O}) > exp{ — (O)}. However, this approach was rejected in the pure 
gauge analyses and we will not pursue it here. For the more interesting gauged SU(2) 
Higgs model the analysis in [26, 27] was quite successful in reproducing the known Monte 
Carlo data but the accuracy obtained fell short of that obtained when using similar 
methods for pure gauge models. This is understandable as the analysis revealed that 
adding the Higgs sector produces new problems (as indeed it does for numerical lattice 
Monte Carlo studies), and the variational ansatz used was relatively simple. 

Given that the Higgs sector causes new problems, it is worth considering pure Higgs 
models on the lattice in the context of OPE methods, effectively the zero gauge coupling 
limit of the gauged Higgs models rather than the infinite Higgs mass limit represented 
by the pure gauge models. However, we have found only a discussion of the real scalar 
model on the lattice using OPE method [12] rather than models with continuous global 
symmetries relevant to realistic gauged Higgs models. 

The purpose of this paper is to present results for the simplest prototype of a Higgs 
sector in a gauge theory, a pair of scalar fields with an 0(2) global symmetry. It is hoped 
that the lessons learned in this case will improve future work on realistic problems. The 
model has not been studied before using the OPE on the lattice and further the Higgs 
sector in this model is different from that used in [26, 27] where only the modulus of the 
Higgs was required. Likewise the variational ansatz we choose here is more complicated 
with more variational parameters than that used in [12] or [26, 27]. We also study different 
schemes to fix the variational parameters. 

Our results show that the OPE can work well for problems with global continuous 
symmetry and large coupling constants. We highlight the role of the symmetry in the 
choice of the variational ansatz and how to use this to reduce the number of necessary 



variational parameters. We also show that minimising the free energy alone is sufficient 
to locate the phase transition. Further using the minimum of the free energy to fix the 
variational parameters for calculations of any quantity produced sensible results. This is 
to be set against our results that show we were unable to find turning points of the form 
(1) for several different types of expectation value. Thus it appears that the principle of 
minimal sensitivity [28] does not work in this model in contrast with its great success in 
pure gauge models [7, 19, 20, 21] and the rigorous results obtained with this approach to 
OPE in other contexts [13, 14, 18]. 

Our starting point is the Euclidean action in four- dimensions for a scalar field theory 
with a global 0(2) symmetry, or equivalently a global U(l) symmetry, put on a four- 
dimensional hyper-cubic space-time lattice to regulate the UV divergences 
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The subscript n runs over all N lattice sites and the n + fi subscript stands for the nearest 
neighbour terms. The aim is to calculate the partition function or generating functional, 
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for large couplings A. We will also be interested in various expectation values such as ((f)) 
and (<p 2 ). 

The OPE method we will use to perform such calculations is the linear (^-expansion 
[7]. A suitable trial action So is added and subtracted from the real action, so that we 
use the action Sg — So + 5(S — So) and 5 = 1, i.e. 

OO r. n 

Z = J2 Z ^(( S o-ST)o (4) 
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where (. . .) is called the statistical average defined as 
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Z is the partition function for the action So. Note that 5 is just a formal counting 
parameter since it must be set to one to make contact with the physical theory. It is 
used to help us truncate the infinite number of Feynman diagrams which represent the 
solutions to QFT, i.e. we will work up to a finite power of 5. It seems that the ultimate 
success of any calculation will rely on the small size of ((So — S) n ) . Thus we might 
guess that, as with the ansatz in any variational approach, S needs to be as good a 
representation of the correct physics as possible. However, this requires us to know the 
full answer before we can make a good guess for Sq. To break this circle, we use our 
physical intuition to make the best guess for the general form of So but at the same time 
we let So depend on some unphysical variational parameters. When we truncate our 
expressions such as (4) at some finite order in delta, our answers will depend on these 
parameters and we will optimize our values for physical quantities by choosing values for 
these unphysical parameters using suitable criteria. Thus we are optimising our choice 
of So using a variational method and it is the optimization which introduces the highly 
non-linear and non-perturbative element of OPE. 



In our case we choose our trial action, So = So((f>, ip; m 2 , A; j^, j^, k+, to be 
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introducing the unphysical variational parameters j^, j^,k+,k-. We do not introduce a 
kfiij) term as we can always make an 0(2) rotation of our fields to eliminate such a term. 
We therefore assume that this has been done so that in general we have used up our 
freedom to make such a field redefinition. Our So is then a reasonable guess for an action 
with two scalar degrees of freedom since the form of So is very similar to S and their ultra 
local parts are identical when — — = and k + = —m 2 . We therefore might hope 
that ((So — S) n )o is likely to be small, and indeed perhaps minimised for values of the 
variational parameters where the ultralocal parts of S and So are equal. Note that our So 
is significantly different from that of [26, 27] as when ^ or either j ^ then this So 
is not 0(2) symmetric. We shall also work with the full complex field (or its equivalent) 
rather than just its modulus as used in the gauged model in [26, 27]. 

However, there is a second criterion for the choice of So, namely it must be of a form 
which allows us to perform some calculations. In this case our trial action So consists of 
ultra local terms. The partition function then factorizes into a product of N integrals 
over the fields at each lattice site. 

Let us first consider the Free energy density, F. Using 5 as a counting parameter and 
define 



where we substituted equation (4) for Z. It is straight forward to equate powers of 5 to 
get a relation between the K^s, the cumulant averages, and the statistical averages of 



The statistical averages required to calculate the K n can then be seen to involve integrals 
over the fields at no more than (n + 1) connected lattice sites. Descriptions of the 
diagrammatic method used to find all the relevant integrals are given elsewhere, for 
instance see [4, 26, 27]. 

For the free energy we choose the unphysical variational parameters, j<j>,j$,k±, such 
that they minimize the free energy (maximise the entropy). For the free energy calculation 
only, it is also the procedure that the principle of minimal sensitivity would demand. After 
fixing the parameters, F depends only on the parameters m 2 and A. 

We also studied the expectation values ((f)), (if>), ((f) 2 ) and (if> 2 ). To do this we again 
use the cumulant expansion, which for the expectation value of an operator O can be 
expressed as 
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Equating powers of 5 produces the same type of relation between statistical averages, 
involving simple integrals over n connected lattice sites, and the L„'s, similar to that for 

A 



K 3 in (8). In practice we expand up to third order, L3 and K3. (0) is then a function 
of the physical parameters m 2 , A and the variational ones, j^, j^, k+, k-, but it does not 
depend on N - the number of lattice sites. 

At this point note that there are alternatives for the optimization step used to set these 
variational parameters in these expectation values. The principle of minimal sensitivity, 
used with great success in many circumstances and in particular in [7, 19, 20, 21], is 
based on (1) where a G {j^, j^, k + , k_}. We applied this method to (0), (ip), (0 2 ) and 
(?/> 2 ). Anticipating the role that the symmetry plays (see discussion below) we also studied 
the 0(2) invariants ((f)) 2 + (tp) 2 and (0 2 ) + (ip 2 )- However we were unable to find any 
suitable solutions for the variational parameters when we tried to minimise the variation 
in the expectation values of the field with respect to the variational parameters. Therefore 
we reverted to the method used in [12], and in calculations of all expectation values we 
choose unphysical variational parameters, j^j^, k±, which minimized the free energy for 
the same set of physical parameters m 2 , A. 

The resulting expressions for each quantity involve a sum of about 100 simple terms, 
each made up of a product of various variational and physical parameters and a few of 
a family of twenty six different integrals. A typical term contribution to the sum which 
makes K 3 in d space-time dimensions is 

I2dj^(m 2 + k + - k_)I 10 I 13 (I 00 y 2 (11) 

where 
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dxdy x a y b exp{ J(j) x + -(k+ + k^)x 2 + j^y + -(k+ - k_)y 2 - \x 2 y 2 }.(12) 

All the integrals are of this form and so they are easily calculated numerically. The 
computational time consuming part was finding the location of the minimum of the 
free energy or of an expectation value to fix the variational parameters. This was done 
numerically in the four- dimensional space of unphysical parameters. Though no numerical 
algorithm can guarantee that it has found the global minimum, in practice it turned out 
that our expressions were amenable to the simplest algorithm. Overall the computations 
took no more than a couple of days on a personal computer. 

It is important to check our results as the expressions are simple but long, so errors 
can easily creep in. One check we used was to set in turn each of the fields to zero. 
The model then reduces to that of a single real field with a 4 interaction and we find 
that we reproduce the results of Wu et al. [12] who used the same approach and an So 
equivalent in this limit. They in turn showed that these results match lattice Monte Carlo 
calculations extremely well. Another check, as we shall discuss below, is to restrict the 
analysis to a subset of one or two variational parameters. If we do that then the results 
agree with the analysis in the full four- dimensional j, k space within numerical errors. 

The results for A = 25 are summarized in the figures. The numerical errors are too 
small to see in the plots but the error inherent in the method can be estimated by the 
difference between the results for the first and third orders. First let us look at the results 
for the free energy as a function of the variational parameters. It is sufficient to study the 
free energy with = and A;_ = as we will show below, so in figure 1 we show results 
for varying and k + for A = 25.0 and m 2 = —15.0. First note that we find no minima or 
any turning points in the even orders. However, figure 1 shows that this is due solely to 
the behaviour of the k + variational parameter which is not always included in previous 
studies.. In the \j\ variable alone, one does seem to find clear minima as shown in figure 

r. 




Figure 1: Plots of the Free Energy from order zero (FO) to three (F3) against j = j<p and 
k = k + for A = 25, m 2 = —15, = and = 0. Contours are evenly spaced in the 
interval -2.2 to -2.55, and show that only the odd orders have minima which get shallower 
at higher orders. 
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Figure 2: Plots of the Free Energy from order zero (FO) to three (F3) against j = for 
A = 25, k + = —m 2 = 15, = k_ = 0. 

2. On the other hand, in the full variational parameter space there are clear minima (and 
not just turning points demanded by the PMS condition (1)) in the first and third order, 
e.g. for A = 25, m 2 = —15.0 these are at k + —0.98m 2 and \j\ pa 1.77. This difference 
in behaviour of even and odd orders is common in OPE methods e.g. appears in [26, 27]. 
Further, the minima in the first and third orders do not move much from the first to 
third order but the minima get shallower at the higher order, i.e. the results become less 
sensitive to the precise value used for the variational parameters. This is characteristic 
of good PMS points. 

Let us briefly consider the principle of fastest convergence of [8] in this context. This 
suggests that one should fix the variational parameters by minimising the difference be- 
tween different higher and lower order terms in the ^-expansion in (7). It is immediately 
obvious that this principle is of no use when we have a function of more than one varia- 
tional parameter, as it provides only too few equations to constrain the several variational 
parameters. However, we can arbitrarily reduce the problem to a one- dimensional one by 
removing the k variational parameters, i.e. set k + = — m 2 and A;_ = so that no change 
in the quadratic part of the action is allowed. These values are also close to the values 
actually found to minimise the free energy. We then set = by symmetry (see below) 
leaving just one free parameter, the source j<p. We are then working with an So similar to 
that used in some of the literature. The results are plotted in figure 2 for example values 
of A = 25 and m 2 = —15 where we are in the spontaneous symmetry breaking regime. 
They show that even in this one-dimensional subspace, the principle of fastest conver- 
gence is still not working well for this Higgs model. The first and third order curves do 
not intersect but the second and third orders do, contradicting normal linear (^-expansion 
behaviour where one compares every second order. However the second and third orders 
do intersect but at a value of much bigger than the location of the minimum. In any 
case the difference between the two orders changes much more slowly around the j = 
intersection 1 suggesting that j = should be chosen by this criterion. It is therefore 

1 Diagramatically, since there are no odd vertices when j — 0, odd orders give zero contributions at 
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difficult to see any reasonable way in which the principle of fastest convergence of [8] can 
be used to fix the variational parameters. 

We do note that the plot in figure 2 does show that the difference between the the 
second and third orders is still small in the region where the Free energy is minimised 
(which is near k + = —m 2 and k- — 0). Thus using the minimum free energy principle 
does give a series with higher order terms giving successively smaller contributions, in 
the spirit of the principle of fastest convergence. 

Now let us analyse the expectation values. One can just use the values for the varia- 
tional parameters which minimise the Free energy, i.e. call upon some sort of maximum 
entropy principle. On the other hand the very successful PMS approach [28] would de- 
mand that we look for turning points in the expectation values with respect to the vari- 
ational parameters (1) and use these values for the particular expectation value needed. 
While that was successful in gauged models [7, 19, 20, 21, 26, 27], we have not been 
able to find any reasonable turning points in the expectation values of the fields or fields 
squared, or even to 0(2) invariant combinations of expectation values. For instance, see 
figure 3. 2 




Figure 3: Contour plot of (0 2 ) + (ip 2 ) against j = and k + at A = 25, m 2 = —15, 
= = 0. 

Thus a major result of this paper is the realisation that the previously successful PMS 
condition does not work for this pure Higgs model. However a minimum Free energy 
criterion does seem to give well behaved results, i.e. weak dependence on the variational 

2 There are some turning points but they occur in regions where the values are changing rapidly. 
There was one shallow saddle point in the (<fi) 2 + (ip) 2 plot at A = 25.0 and m 2 = —15, near fc_ = 25, 
\j\ = 2. However no other expectation value at any order showed any good behaviour in this region and 
this suggests that this was a feature of increasing ripples in the function as k + increased. 



parameters near the chosen solution, small corrections in the variational parameters and 
free energy values when moving from first to third orders. There is no reason to believe 
that this would not be the case in gauged models, at least for weak gauge coupling. 

Now let us try to use the variational values obtained from the minimal Free energy 
principle to study the physics in this model for the large coupling A = 25. The most 
prominent feature of the complex 4 theory is the phase transition. Figure 4 shows a plot 
of the quantity ((f) 2 + (if)) 2 vs. m 2 . We can clearly see that the curves are continuous 




Figure 4: This plot shows the first and third order result of ((f)) 2 + (if)) 2 vs. m 2 for A = 25. 
The phase transition occurs around m 2 « —13.4 for third order and m 2 « —11.8 for first 
order. 



but that there is a sudden change in slope at m 2 ~ —13.4 for the third order expansion 
(and at m 2 ~ —11.8 for the first order expansion). This indicates a second order phase 
transition. Since we are working on a lattice and have not taken the continuum limit, the 
mass parameter m is not a physical mass so the fact that the phase transition does not 
occur at m 2 = (as it would classically) is not surprising. In any case we are working 
a long way from the perturbative regime. The fluctuations about the minimum can be 
estimated from the difference {(<f> 2 ) — ((f) 2 ) + ((ip 2 ) — (4>) 2 ) an d this is shown in figure 5. 
This peaks at the phase transition as expected. 

Another point to note is the way that non-zero results for ((f)) and (if>) can be found in 
the OPE. This is a key difference between this method and lattice Monte Carlo methods. 
The OPE will always find one and only one vacuum solution even if there is more than 
one solution. It does this by allowing the sources j to take non-zero values in which case 
the classical potential is tilted favouring one vacuum solution. This contrasts with lattice 
Monte Carlo methods which sample all possible vacuum states with equal likelihood. This 
means OPE methods are ideal for studying non-trivial field configurations, perhaps with 
defects present, by using classical sources to manipulate the effective quantum potential. 

The last plot, figure 6, clearly shows the role of the 0(2) symmetry in this method. 
The plot of (if>) vs. (<f)) for m 2 = —15.0 (inner curve) and m 2 = —18.0 (outer curve) both 
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Figure 6: Plot of (ip) vs. ((f)) for m 2 = —15.0 and m 2 = —18.0, for various initial values 
for the variational parameters. A = 25.0. We see that there are many degenerate vacuum 
states for a given m 2 in the broken symmetry regime. 
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show circles, i.e. the quantity (<f)) 2 + (ip) 2 is a constant. 

It is important to note the form of the solutions found for the unphysical parameters 
j(f>,jip,k+ and In the unbroken phase, = = 0. In the broken phase while 
\j\ 2 = 3 6 + i| i s a non-zero constant for a given set of physical parameters m 2 and A, 
with a different solution for the j's corresponding to each point on the circles in figure 
6. This immediately suggests that \j\ acts as an order parameter. Thus a calculation of 
just the free energy, and by implication the variational parameters, rather than of any 
expectation value can be used find the location of the transition, as shown in figure (7). 
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Figure 7: Plot of \j\ vs. m 2 for A = 25.0. It shows that this variational parameter acts 
as an order parameter. 

The solution for k + was equal to m 2 in the unbroken regime and within 10% of m 2 
otherwise, as shown in figure 8. Thus (K + + m 2 ) is acting as another order parameter. 
This shows that its numerical value, chosen by minimising the free energy, is keeping S 
and 5*o as similar to each other as possible. However, it is not so much the numerical 
value of k + but the behaviour of functions as we vary k + which is important. As remarked 
above, there is often no minimum in the k + variable when there are minima in the other 
variational parameters. 

Finally we note that that k_ = always minimised the free energy, within numerical 
accuracy. This latter result is of particular significance with regards the symmetry of our 
ansatz So, and we will now turn to consider this aspect. 

To explain these results, in particular figure 6, we must consider the symmetry of the 
problem. S is invariant under 0(2) transformations such as 

^(- c :% *=(*)• (i3) 

Our So of (6) is not 0(2) invariant if ^ or if the j parameters are non-zero. However 
it is easy to see that 

So(0; m 2 , A; j, k + , k. = 0) = S $ = U$; m 2 , A; f = Uj, k + , k_ = 0) (14) 




Figure 8: Plot of (—1 — k + /m 2 ) which minimise the third order free energy vs. m 2 for 
A = 25.0. 

where j = (j'^iv)- That is, there is a set of field independent transformations for j^jj^ 
and k + which leave So invariant under the 0(2) field transformation. Note that we must 
have = for this to be true, so let us for the moment assume that this value minimises 
the free energy 3 . 

Using this 0(2) transformation we note that the integration measure is also 0(2) 
invariant. It is then easy to show that at whatever order in S we truncate our free energy 

expression, if the set {j, k + , = 0} minimises the free energy then the set {j', k + , k^. = 0} 
also minimises the free energy. It then follows that a calculation of any 0(2) invariant 
physical object will not depend on which solution for the variational parameters we use. 
In particular the free energy, ((f) 2 ) + (ip 2 ), and ((f)) 2 + (if)) 2 will all give results which are 
independent of the solution for the variational parameters. For objects like (<fi) which 
are not 0(2) invariant, we can still relate one solution found for it, with a particular 
set of variational parameters, to another solution for ((f)) and another set of variational 
parameters. The relationship is precisely the simple 0(2) transformations discussed here. 

Plotting out final values for the variational parameters and for the different 
initial conditions used in producing figure (6) we found that and also form circles 
in a j^-j^ plane (for fixed m 2 and A). Actually, there is a one-to-one correspondence 
between pairs of j</,, and ((f)), (if)) - two points separated by a certain angle on the 
j circle correspond to two points separated by the same angle on the expectation value 
circle; so to get one of the points from the circles in figure 6, we fix to a certain 
value, and minimize F in the remaining 3-dimensional parameter space. Once the point 

3 Remember in our choice of So in (6) we used the 0(2) field rotation to eliminate a possible k<jytp 
term. However, with L = we regain the ability to exploit this freedom. 
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is found, we fix to a different value, and minimize again, thereby finding a different 
point. This is the explicit display that 0(2) symmetry transforms both j's and 0's with 
the same rotation matrix. The fact that we found k + was constant for all these solutions 
and that k_ = confirms the role of the symmetry transformations (13) and (14). 

The presence of the variational parameter which does not respect the 0(2) sym- 
metry even if we allow field independent transformations, is a crucial distinction from the 
work of [26, 27] and [12]. It means that we have allowed our model to break the bounds of 
the symmetry of S if the dynamics so choose. Thus one of the most important numerical 
results we have is that = to numerical accuracy (at least six significant figures). 
Once this is known, only then do the symmetry arguments of the previous paragraph 
explain the circles of figure 6. 

This explicit display of the role of the symmetry in the solutions suggests an interesting 
check, namely to impose the symmetry on the Sq. Thus, we ran the whole process again 
holding A;_ = 0, i.e. minimizing in a 3-dimensional parameter space. We recovered the 
circle of j's (at fixed m 2 ). We can then take it a stage further as in addition to holding 
the k- = 0, we can also fix one of the j's to zero, and thus minimize in a 2-dimensional 
parameter space. In the un-broken symmetry regime, the j's are zero anyway, while in the 
broken symmetry phase there is an infinite set of equally valid solutions for the variational 
parameters, so holding one of the j's at zero does not imply any loss of generality of the 
results. 

In principle, one of the big advantages of the OPE is that one can systematically study 
higher orders through a straight forward extension of the method used here. However, it 
is likely that no solution for the variational parameters will be found at the next order, 
just as none was found at second order. As mentioned earlier, this behaviour is common 
in an OPE. Thus we would need to calculate K§ or L5.. 

To summarize, the OPE method used here works well in identifying the phase transi- 
tion when there is a continuous symmetry, provided we fix the variational parameters by 
minimizing the free energy and we do not use the previously successful principle of min- 
imal sensitivity. The explanation for the failure of the latter approach in this case after 
its earlier successes seems to be our use of a variational parameters which are quadratic 
in the Higgs fields. Another clear message coming from this work is that the optimal 
solutions choose a trial action So which has the same symmetry as the full action S. This 
fact can be exploited to dramatically reduce the number of variational parameters used, 
which in turn will simplify and accelerate the analysis of more complicated models. Next, 
the source variational parameters, j, and the combination (k+ + m 2 ), act as order pa- 
rameters. Thus one can identify the transition point from free energy calculations alone 
which is a further simplification when trying to find the phase diagram. We also note 
that this method can be used to study just one of many equivalent vacuum solutions 
which distinguishes it from lattice Monte Carlo. On the other hand, as shown here with 
just three orders, the OPE has a practical and systematic way of improving its accu- 
racy, unlike other analytic non-perturbative methods. Thus we have demonstrated that 
OPE offers a practical route to the study of models with Higgs sectors and continuous 
symmetries in non-perturbative problems. 
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Figure 9: This plot shows the first and third order result of ((f) 2 ) + (ip 2 ) vs. m 2 for A = 25. 
It shows the phase transition, but the kink is not as clearly visible as in other plots. 
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